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TIME-DEPENDENT SLOW-TIME-SCALE THEORY 
OF FREE-RUNNING AND PHASE-LOCKED 
GYROTRON OSCILLATORS 


I. INTRODUCTION 

There is currently considerable interest in the development 
of high power phase-locked gyrotron oscillators. These devices 
have the potential to combine the high efficiency and power 
associated with oscillators with the coherence and phase control 
properties associated with amplifiers. Although previous 
theoretical work on gyrotron ope ration!1-8] has been successfully 
applied to the development of cw devices for heating of tokamak 
plasmas, the investigation of phase-locked gyrotron operation, 
particularly involving pulse-line accelerator driven devices, 
presents a new challenge to the theorist. These devices, based 
on intense relativistic electron beams, are characterized by 
short pulselengths and voltage risetime and ripple effects which 
can prevent the establishment of a steady-state equilibrium. A 
consideration of non-stationary operation is also involved in the 
investigation of conditions for phase-locked operation. Time- 
dependent effects can be studied using a particle-in-cell 
simulation code of the type developed by Lin and co-workers[9]. 

In this work an alternative slow-time-scale approach is used to 
obtain a relatively simple, economical analysis of the time- 
dependent, nonlinear dynamics of pulsed gyrotrons. 

The time-dependent theory of gyrotrons has been considered 
by Nusinovich and co-workers[10,11], mainly in the context of 
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multimode operation and mode stability. A time multimode theory 
of quasi-optical gyrotrons has been developed by Bondeson et 
al.[12]. Early work on mode selection and phase-locking of 
vacuum tube oscillators was carried out by Van der Pol[13] and 
Adler[14]. An analytical theory of the conditions for phase¬ 
locking gyrotrons has been presented by Manheimer[15,16]. The 
emphasis in the present work is on the modeling of transient 
effects in gyrotrons whose pulse length is only a few times 
longer than the characteristic risetime of cavity fields. The 
approach is developed for the case of several interacting modes 
with closely spaced frequencies and interacting at the same 
harmonic of the cyclotron frequency. Calculations are presented 
for single mode operation of relativistic gyrotron with 
parameters similar to an experimental device under test at the 
Naval Research Laboratory(17]. Both free-running and phase- 
locked operation are considered. The phase-locking signal is 
introduced via a prebunched beam. The effect of a voltage ripple 
is investigated. The possibility of phase-locked operation in 
the presence of voltage ripple - an important issue for proposed 
NRL experiments - is demonstrated. Calculations for multimode 


operation will be the subject of future work. 

The time theory of gyrotron oscillators derived in this work 
is based on slow-time-scale equations for the electron motion 
similar to those used in steady-state models[7,18,19]. Slow- 







time-scale equations for the cavity rf field amplitude are 
obtained by expressing the time-dependant behavior relative t 


reference frequency a> 0 which is close to the operating freque 
«. The next section of this paper contains an outline of the 
theoretical approach. The third section contains the results 


calculations and the last section presents some conclusions d 
from this research. 









II. THEORY 


Consider a gyrotron with a cylindrical resonator and a th 
annular beam. The electrons follow helical trajectories in thi 
applied axial magnetic field about guiding centers located at . 
radius R 0 from the symmetry axis. The electron beam interacts 
with one or more resonator modes, which are assumed to be TE 
modes with closely spaced oscillation frequencies. It is 
convenient to look at time effects which remain after a referei 
frequency oo a has been factored out. These effects are 
characterized by time scales which are much longer than the wa< 
period and are incorporated in time mode amplitudes A n (t). Thi 
total transverse electric field is expressed in the form 


E t = l A n (t) h n ( z )e n ( r, 0? z )exp(-iw 0 t) (1) 

n = 1 

where e n is a transverse mode vector function defined in Ref. 
and h n (z) is an axial profile function which satisfies the 
Helmholtz equation; 









The function P n »h n e n is an eigenfunction of the empty resonator 
in the single transverse mode approximation with eigenfrequency 
<*> no . Thus the present formulation neglects beam loading effects 
on the mode axial profile. N is the number of interacting modes 
The slowly varying mode amplitudes have the form: 

A n (t) - f n (t) e -1 ^" (t) (3) 

The phase derivative ^“d^/dt corresponds to the frequency 
difference between the mode oscillation frequency w n and co 0 . 

The transverse electric field satisfies the wave equation: 

. 3 2 E t 3 J t 

V 2 E t - -- =* fj 0 - ( 4 ) 

c 2 3t2 3t 

where J t is the transverse ac current density, c is the speed of 
light and fj Q is the permeability of free space. MKS units are 
used throughout except as noted. Substituting Eq.(l) into 
Eq.(4), introducing a mode Q factor, Q n , to account for 
diffraction and ohmic losses, and noting that dA n /dt<<w 0 A n one 













where co no is the empty resonator eigenfrequency for the nth mode. 
Multiplying Eq.(5) on both sides by (h n e n )*, integrating over the 
cavity volume and applying the orthogonality properties of the 
transverse vector functions!18] leads to the mode excitation 
equation: 


[-♦ — 

L dt 2Q 


i ( w 0 - w n ) A n (t) 


ic 2 u a r 9J t 

- dzda h* el •— e la) 


where N n is the axial mode normalization integral: 


N n - J L dz |h n (z) | 2 


To evaluate the right hand side of Eq.(6) the rf current density 
is written in the approximate form(7J: 


J t =J exp( -ioo 0 t) 


( 8 ) 


p 2 n £ 

= d(w 0 t)J t e 

o J 0 


Introducing the normalized frequency shift & n = ( w 0 -<*> n 0 )/(*> 0 and 
time t=w 0 t and substituting Eq.(8) into Eq.(6) leads to: 


d 1 , c* //, - 

— + - - i S n ( t ) = - —- Jdzda h* e‘ • J (j) e* 1 


dT 2 Q n 
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To calculate the ac current density, the interaction with the 
electron beam is treated in the single particle approximation. A 
considerable simplification of the general time-dependent problem 
results if one uses the fact that the characteristic rise-time of 
fields in the resonator is much longer than the electron transit 
time in the cavity. In this case one can use a quasi-steady- 
state approximation in which the electron trajectories are 
calculated for rf fields with fixed amplitude, .:), and 
linearized phase, '( t 0 )( t-t 0 ). The slow-tiroe-scale 

nonlinear electron equations of motion for a particular harmonic 
are readily deduced from previous steady-state analyses[7,18] and 
are given by: 


du t 

dz 




dhn ] c -i[ A+^-Qn.-s)*,. ] 
dz ' 


(11a) 


dA sy M sJ s (k n r L ) u z dh n wj( 0 u£ . 

-3- - - - Z F n —3—--Re I fh n + i —-3- ——h n j 

dz u z u t n = 1 k n r L P y« 0 dz s$2&> 0 r 

e -i[ A+* n -(m„-s)3 0 ]| + 5 r x _ Jfla ) U1 

J (0 o r 

du z u t n r dh n 

— = —— Z F n J'(ii n r L )Re i — 
dz u z w 0 n = 1 L dz 


e ~i [ A+\j/ n - (m n -s ) Z 0 ] 


(11c) 













where u t »yv t /c is the normalized transverse momentum amplitude, 
u z -yv,/c is the normalized axial momentum, A is a slowly varying 
transverse momentum phase, s is the harmonic number, m n is the 
mode azimuthal index, y is the relativistic mass ratio, k n is the 
mode transverse wavenumber, r L is the Larmor radius of the orbit, 
J s (J' s ) is (the derivative of) a regular Bessel function, 5 is 
the relativistic cyclotron frequency correspondong to the input 
beam energy, S 0 is the orbit guiding center azimuthal angle, and 
F n is the normalized rf field amplitude: 



Quantities with a " have been normalized according to: 
z-z/r w0 , r L -r L /r W0 , 5-$?r MO /c, u> 0 r w0 /c, and E n -k n r wo . R„ 

denotes the orbit guiding center radius, e is the electron 
charge, m # is the electron mass, x' n is a zero of j; for the mode 
n, and r wo is an arbitrary normalization length usually chosen to 
be a representative cavity wall radius. 

The ac current density is obtained by integrating Eqs.(ll) 
for an appropriate set of initial conditions at the cavity input 
at z 0 . For a cold, phase-mixed electron beam: u t (z 0 )=u t0 , 


u z (z 0 )= u j. 0 , and A(z 0 )=A 0 is uniformly distributed in the interval 
[0,2n]. The case of an electron beam prebunched by passage 









through a bunching cavity driven 
by an initial phase angle distrib' 

A o = 0 o + q B sin0 o - 

where 0 O is uniformly distributed 
the bunching parameter, and * 0 is 
annular beam, Eqs.(ll) are averag 







l e ko 


B ( k„ R 0 ) 


m Q cw 0 n( l-m 2 /*^ 2 ) Jjj; )N n 


and < > a denotes an average with respect to the variable a. The 
time behavior of a free-running or driven gyrotron oscillator is 
obtained by solving Eqs.(15) together with (11) subject to 
appropriate initial conditions as discussed above. When the 
phase \|/ n of a driven oscillator (via beam prebunching) converges 
to a constant value, the oscillator is said to be phase-locked. 
The locking frequency is w 0 in the present formulation. When the 
interaction occurs with only a single circularly polarized mode, 
the ac current density averaged with respect to A e is independent 
of the guiding center azimuthal position so that the average with 
to S 0 may be omitted. 
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III. RESULTS 


Calculations have been carried for a high voltage 35 GHz 
gyrotron configuration similar to the NRL experiment recently 
reported by Gold et al.[17]. The peak voltage and current of the 
annular beam are taken to be 650 kv and 1.5 kA, respectively. 

The beam guiding center radius is R 0 -1.16 cm, the cavity radius 
is 1.6 cm, and the longitudinal profile of the cavity fields is 
assumed to be sinusoidal with length 4 cm. The operating mode is 
the TE 621 circular mode with polarization counter-rotating to the 
beam rotation. The beam pitch ratio asv t /v z »l. The cavity Q 
factor Q-250. The applied magnetic field is 2.4 Tesla. Spreads 
in beam guiding center and pitch ratio are neglected as are 
space-charge effects. In these calculations the cavity is seeded 
with low amplitude rf fields (1 kV/cm at the beam in the present 
case). The results of calculations for free-running and phase- 
locked oscillator operation are presented below. The effect of a 


voltage ripple is shown. 

The first set of calculations correspond to an idealized 50 
nsec voltage pulse with a 3 nsec rise from 75% of peak to peak 
voltage and no ripple. The beam current and pitch ratio are 
assumed to vary linearly with voltage and as V 1 - 33 , respectively, 
during the voltage rise up to the maximum values given above. 
Results for the free-running oscillator are shown in Figures 1 
(a)-l(c). These Figures show, respectively, the electronic 
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efficiency, operating frequency, and rf phase of the oscillator 
plotted as a function of time. Figure 1(b) plots the shift in 
operating frequency oo relative to w 0 which is taken to be the the 
cold cavity eigenfrequency (which is close to 35 GHz). These 
plots indicate that steady-state operation is achieved after 
about 8 nsec. The steady-state efficiency is 17.7%, 
corresponding to an output power of 173 MW, and the steady-state 
frequency shift 8 rR -(oi rR -<»> 0 )/w 0 due to beam loading is 5.5xl0 -4 
which corresponds to 19.2 MHz. This frequency shift leads to a 
linearly increasing rf phase with time as shown in Figure 1(c). 

Results for an oscillator driven via a prebunched beam are 
given in Figures 2(a)-2(c). The beam is modulated at a frequency 
w 0 and the modulation amplitude is modeled by the bunching 
parameter q B =0.5. An estimate for the locking signal power 
required to obtain a beam with q B -0.5 is given below. In this 
calculation the difference between the normalized locking signal 
frequency and the normalized free-running oscillator (FRO) 
frequency is 8 L = ( w 0 -<*) F R )/w 0 =1.05x10" 3 or Af-35 MHz. Figures 2(b) 
and 2(c) show that phase-locked operation is obtained after 30 
nsec into the voltage pulse. The phase-locked condition is 
indicated by a frequency shift <*>-w o =0 and constant rf phase. As 
shown in Figure 2(a), in this case the oscillator efficiency is 
reduced slightly to 16% while phase-locked. 













If the difference between the locking signal frequency and 
the FRO frequency Af is too large, phase-locking will not occur. 
Instead a periodic perturbation of the oscillator behavior is 
observed as shown in Figures 3(a)-3(c). The locking signal 
detuning in this case is 8 L »1.55xlO~ 3 or Af-54 MHz. The beat 
frequency of the perturbation is about 28 MHz. The non-phase- 
locked behavior of the rf phase is clearly shown by Fig. 3(c). 

The case of a FRO with a voltage ripple is shown in Figures 
4(a)-4(c). The voltage ripple amplitude is 2% and the ripple 
period is 10 nsec. Other parameters are the same as for the FRO 
case discussed above. Figures 4(a) and 4(b) show that this 
voltage ripple causes a 3.4% amplitude ripple in the efficiency 
and a 0.06% amplitude ripple (21 MHz) in the frequency. Figures 
5(a)-5(b) show the results of driving this oscillator 
configuration with a prebunched beam with q„=0.5 and modulation 
frequency shift -5.5x10“ 4 (19 MHz). Figures 5(b) and 5(c) sho 

that phase-locked operation is achieved in the sense that after 
about 30 nsec the oscillation frequency averaged over the ripple 
period equals the drive frequency and the averaged rf phase angl 
is a constant. As in the FRO, the voltage ripple leads to a 
0.06% frequency ripple in the PLO. The frequency ripple causes 
corresponding 13° amplitude ripple in the phase angle. Figures 
6(a)-6(c) represent the case of an unphase-locked, driven 
oscillator with a 2% voltage ripple. The drive signal frequency 










shift is $ l -2.5x10“ 3 (88 MHz). The beat frequency for the 
perturbed oscillator is about 34 MHz. The failure to achieve 
phase-locked operation is clearly shown by the phase angle 
plotted in Figure 6(c). 

To estimate the input power needed to obtain a beam bunching 
parameter q B «0.5 for circularly polarized radiation and a 
gaussian rf profile in the bunching cavity, one can apply the 
small signal gyroklystron theory developed by Tran et al [20]. 

The required input power in Watts is: 

4itc 3 m 3 <x' 2 -m 2 ) J 2 (x' n ) rl0? o 8 zo q| e 2x2 

P in (W) - -2- 

24 2fj Q e 2 J 2 (k n R 0 ) fj x (4 3/t/j+^/j) Q x 

(17) 

where -n(3 2 t 0 L t /& z 0 X is the normalized interaction length of the 

bunching cavity, is the normalized drift length, Q x is the 
bunching cavity Q factor, and x=( w 0 -S)L l /4v J!0 the bunching cavity 
magnetic field detuning parameter. The other parameters are 
defined above. Typical choices for // x and y A are 2 and 3, 
respectively. Assuming the same mode, magnetic field and beam 
parameters in the bunching and power cavities, and Q 1 =200, one 
obtains P in =>530 kw as the power required to phase-lock a gyrotron 


oscillator with an output power of over 160 MW for an effective 
gain of over 24 dB. 












IV. CONCLUSIONS 


A time slow-time-scale theory and a computer code for 
numerical calculations have been developed to investigate pulsed 
free-running and phase-locked gyrotrons. This code is expected 
to be a valuable resource for the design and interpretation of 
planned phase-locked gyrotron experiments. The calculations 
presented here show the feasibility of observing phase-locked 
operation with 50 nsec pulses. Phase-locking over a significant 
bandwidth is shown to be possible with a locking power 24 dB less 
than the oscillator output power for a single prebunching cavity. 
The effect of voltage ripple has been investigated and shown not 
to prevent phase-locked operation, however, applications 
requirements for rf phase control may place significant limits on 
allowable voltage ripple. 
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